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ABSTRACT 

In this paper we study the stellar-mass dependence of galaxy clustering in the 6dF 
Galaxy Survey. The near-infrared selection of 6dFGS allows more reliable stellar mass 
estimates compared to optical bands used in other galaxy surveys. Using the Halo 
Occupation Distribution (HOD) model, we investigate the trend of dark matter halo 
mass and satellite fraction with stellar mass by measuring the projected correlation 
function, w p (r p ). We find that the typical halo mass (Mi) as well as the satellite power 
law index (a) increase with stellar mass. This indicates, (1) that galaxies with higher 
stellar mass sit in more massive dark matter halos and (2) that these more massive 
dark matter halos accumulate satellites faster with growing mass compared to halos 
occupied by low stellar mass galaxies. Furthermore we find a relation between Mi and 
the minimum dark matter halo mass (M m i n ) of Mi ea 22M m ; n , in agreement with 
similar findings for SDSS galaxies. The satellite fraction of 6dFGS galaxies declines 
with increasing stellar mass from 21% at M st ciiar = 2.6 x 10 10 h~ 2 M Q to 12% at 
Mgtciiar = 5.4 x 10 10 h~ 2 Mq indicating that high stellar mass galaxies are more likely 
to be central galaxies. We compare our results to two different semi-analytic models 
derived from the Millennium Simulation, finding some disagreement. Our results can 
be used for placing new constraints on semi-analytic models in the future, particularly 
the behaviour of luminous red satellites. Finally we compare our results to studies of 
halo occupation using galaxy-galaxy weak lensing. We find good overall agreement, 
representing a valuable crosscheck for these two different tools of studying the matter 
distribution in the Universe. 

Key words: galaxy formation, large-scale structure of Universe, surveys, galaxies: 
statistics, galaxies: halos 



1 INTRODUCTION 

The first statistical studies of galaxy cluster 



ing jTotsuji fc Kiharalll969l ; |Peebledfl~973l ; lHauser fc Peebles! 



Il973l . [l9r 



Peebles! 1974 ) found that the galaxy correlation 
function behaves like a power law, which is difficult t o 
explain from first principles (jBerlind fc Weinberg! l200ff ). 
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More recent studies, h owever, found devia tions from a 
power law. For example IZehavi et al.l (|2005al ) showed that 
the projected correlation function w p (r p ) of SDSS galaxies 
exhibits a statistically significant departure from a power 
law. They also showed that a 3-para meter Halo Occupa- 



tion D istribution (HOD ) model (e.g.. iJing. Mo fc B orner 



19981 : iMa fc Fry! hood; iPeacock fc Smith! l200d: ISeliak 



200d: IScoccimarro et all l200ll : iBerlind fc Weinberg! TSoOll 



Coorav fc Shetrj|2002l ) together with a ACDM background 
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cosmology, can account for this departure, reproducing the 
observed w p (r p ). 

Within the halo model the transition from the 1-halo 
term to the 2-halo term causes a "dip" in the correlation 
function at around 1 — Sh^ 1 Mpc, corresponding to the ex- 
ponential cutoff in the halo mass function. In case of a 
smooth transition between the one- and two-halo terms, this 
can mimic a power-law correlation function. Studies with 
Luminous Red Galaxies (LRGs) found that the deviation 
fro m a power-law is large r f or highly clustered bright galax 



trom a power-law is larger lor highly clustered bright galax- 
ies (IZehavi et al.l l2005aHbl; iBlake. Collister fe LahavlBooi : 



IZheng et al.l |2009|; IZehavi et al.l |2010T). an d at high red 



shift (|Conrov. Wechsler fc Kravtsov 12006 ), wh ich agrees 



with theoretical predictions (jWatson et al 



201ll) . 



While galaxy clustering is difficult to predict, dark mat- 
ter clustering is dominated by gravity and can be predicted 
for a given cosmology using N-body simulations. Using mod- 
els for how galaxies populate dark matter halos, which are 
usually motivated by N-body simulations, we can directly 
link galaxy clustering and matter clustering. This can be 
modelled in terms of the probability distribution p(N\M) 
that a halo of virial mass M contains N galaxies of a given 
type. On strongly non-linear scales the dark matter distribu- 
tion is given by the actual density distribution of the virial- 
ized halos, while on large and close to linear scales the dark 
matter distribution can be predicted from linear perturba- 
tion theory. 

HOD modelling has been applied to galaxy clus- 
tering data from the 2-degree Field Galaxy Redshift 
Surv ey (2dFGRS) dPorciani. Magliocchetti fc Norberel 
120071 ; iTinker et all 120061) and the Sloan Dig ital Sky 
Survey (SPSS) Ivan den Bosch, Yang fc Moi 12003 



Magliocchetti fc Porcianil |2003| ; IZehavi et all l2005allr 



Tinker et all 120051 ; lYang et all [20051 . 120081 ; IZehavi et ail 
2010h . More recently it also became possib le to model 



the c lustering of high-z galaxies using VV DS (jAbbas et all 
|2010| ). Bootes (|Brown et al.ll2008l l DEEP2 (|Coil et allbOOrj ) 
and Ly man-break galax ies at high redshift in the GOODS 
survey l|Lee et all 2006). Such studies revealed, that the 
minimum mass, M m j n for a halo to host a central galaxy 
more luminous than some threshold, L is proportional to L 
at low luminosities, but steepens above L*. Massive halos 
have red central galaxies with predominantly red satellites, 
while the fraction of blue central gala xies increases with de - 
creasing host halo mass. Furthermore IZehavi et all (|2005bl ) 
found that there is a scaling relation between the minimum 
mass of the host halos, M m i n and the mass scale, Mi of 
halos that on average host one satellite galaxy in addition 
to the central gala xy, Mi » 23M m j n . Using a different HOD 



parameterization, IZheng. Coil fc Zehavil d2007l) found the 
relatio n to be Mi « 18M m i n , very similar to TZehavi et all 
(|2010l ) who found Mi « 17 M min . 

The 6dF Galaxy Survey is one of the biggest galaxy sur- 
veys available today with a sky coverage of 42% and an aver- 
age redshift of z — 0.05. The survey includes a bout 125 000 
redsh i fts se l ected in the J, H, K, bj, rF-bands (|jones et all 
120041 . I2005I . l2009l l The near infrared selection, the high 
completeness and the wide sky coverage make 6dFGS one 
of the best surveys in the local Universe to study galaxy 
formation. This dataset has been used to study the large 
scale galaxy clustering to measure the Hubble constant us- 
ing Baryon Acoustic Oscillations (|Beutler et all l201ll ). as 



well a s the growth of structure at low redshift jBeutler et all 
|2012|) . While these previous studies used the Tf-band se- 
lected sample, in this analysis we use the J-band. The J- 
band allows the most reliable stellar mass estimate of the 
five bands available in 6dFGS, because of its lower back- 
ground noise. Together with the bj — rp c olour we can de- 
rive s tellar masses using the technique of iBelT fc De Jongl 
(2001). which leads to a dataset of 76 833 galaxies in to- 
tal. The photometric near-infrared selection from 2MASS 
makes the stellar mass estimates in 6dFGS more reli- 
able than stellar mass estimates in other large galaxy 
surveys which rely on optic a l bands (iDrorv et all 
Kannappan fc Gawisei 20071; Longhetti fc Saraccol 



2004 



2009: 



Grillo et al.ll200Sl ; iGallazzi fe Bellll2009h 

Numerical N-body simulations are usually restricted to 
dark matter only. To understand galaxy formation, bary- 
onic effects such as feedback and gas cooling, have to be 
included. Such simulations face severe theoretical and nu- 
merical challenges. Semi-analytic models build upon pre- 
calculated dark matter merger trees from cosmological sim- 
ulations and include simplified, physically and observation- 
ally motivated, analytic recipes for different baryonic effects. 
Semi-analytic models have been shown to successfully repro- 
duce observed statistical properties of g alaxies over a large 
range of galaxy masses and redshifts (e.g . Croton et al.ll2006l; 
Bower et all 12006); |p~e Lucia fc BlaizotJ 120071 ; iBertone et all 
20071 ; iFont et all boost iGuo et all l2010t ) and allow a level 
of understanding unavailable in N-body simulations. The 
underlying models are necessarily simplified and often use 
a large number of free parameters to fit different observa- 
tions simultaneously. In this paper, we derive 6dFGS mock 
surveys fro m semi-analytic mod els based on the Millennium 
Simulation (|Springel et alll2005l ) and compare the properties 
of these surveys with measurements in 6dFGS. Our results 
can be used to improve upon these semi-analytic models and 
further our understanding of baryonic feedback processes on 
gala xy clustering. 

iMandelbaum et alj (|2006l ) studied halo occupation as a 
function of stellar mass and galaxy type using galaxy-galaxy 
weak lensing in SDSS. They found that for a given stellar 
mass, the halo mass is independent of morphology below 
Msteiiar = 10 11 Mq, indicating that stellar mass is a good 
prox y for halo mass at that ra nge. We compare our results 
with IMandelbaum et all (|2006T ) which represents a valuable 
crosscheck of the HOD analysis using two very different tech- 
niques, weak lensing and galaxy clustering. 

This paper is organised as follows: First we introduce 
the 6dF Galaxy Survey in section [5] together with the tech- 
nique to derive the stellar masses. We also explain how we 
derive our four volume-limited sub-samples in stellar mass 
and redshift, which are then used for the further analysis. 
In section [3] we calculate the projected correlation function, 
w p (r p ), for each sub-sample and use jack-knife re-sampling 
to derive the covariance matrices. In section [3] we fit power 
laws to the projected correlation functions of the four sub- 
samples. In section[5]we introduce the HOD framework and 
in section [B] we apply the HOD model to the data. In sec- 
tion [7] we derive 6dFGS mock samples from two different 
semi-analytic models, which we then compare to our results 
in section [H] together with a general discussion of our find- 
ings. We conclude in section [9] 

Throughout the paper we use r to denote real space sep- 
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arations and s to denote separations in redshift space. Our 
fiducial model to convert redshifts into distances is a flat 
universe with Q m — 0.27, w = — 1 and fi*, = 0. The Hubble 
constant is set to Ho = WOh km s _1 Mpc _1 which sets the 
unit of stellar masses to h~ 2 MQ, while most other masses 
are given in h~ Mq . The HO D model uses cosmolo gical pa- 
rameters following WMAP7 jKomatsu et al.ll201ll ). 



2 THE 6DF GALAXY SURVEY 

The 6dF Galaxy Survey C6dFGS; [jones et al.|[20ol . 120051 . 
120091 ) is a near-infrared selected (J, H, K) redshift survey 
covering 17 000deg 2 of the southern sky. The J, H and K 
surveys avoids a ±10° region around the Galactic Plane 
to minimise Galactic extinction and foreground source con- 
fusion in the Plane. The near-infrared photometric selec- 
tion was based on total magnitudes from the Two-Micron 
All-Sky Survey Ex tended Source Catalog (2MASS XSC; 
IJarrett et al.l |2000| ). The spectroscopic redshifts of 6dFGS 
were obtained with the Six-Degree Field (6dF) multi-object 
spectrograph of the UK Schmidt Telescope (UKST) between 
2001 and 2006. 

The 6dFGS J-selected sample used in this paper con- 
tains (after completeness cuts) 76 833 galaxies selected with 
9.8 J ^ 13.75. We chose the J-band because it has the 
highest signal-to-noise of the three 2MASS bands. While 
there is slightly less extinction in the K-ha.nd compared to 
the J-band, for practical purposes, the J-band has better 
S/N because the night sky background glow is much less 
in the J- than in the if-band. The near infrared selection 
makes 6dFGS very reliable for stellar mass estimates. 

The mean completeness of 6dFGS is 92 percent and 
the median redshift is z — 0.05. Completeness corrections 
are derived by normalising completeness-apparent magni- 
tude functions so that, when integrated over all magnitudes, 
they equal the measured total completeness on a particular 
patch of sky. This pr ocedure is outlined in the luminosity 
function evaluation of Ijones et alj (|2006l ) and a lso in Jones 
et al . , (in prep) . The original survey papers (|jones et al.l 
|2004 120051 . 120091 ') describe in full detail the implementation 
of the survey and its associated online database. 

The clustering in a galaxy survey is estimated relative to 
a random (unclustered) distribution which follows the same 
angular and redshift selection function as the galaxy sample 
itself. We base our random moc k catalogue gener ation on 
the 6dFGS luminosity function (| Jones et al.ll2006l ). where 
we use random numbers to pick volume-weighted redshifts 
and luminosity function-weighted absolute magnitudes. We 
then test whether the redshift-magnitude combination falls 
within the 6dFGS J-band faint and bright apparent magni- 
tude limits (9.8 ^ J ^5 13.75). We assigned a bj-rF colour to 
each random galaxy using the redshift- bj-rF colour relation 
measured in the data and used these to derive stellar masses 
for the random galaxies using the same technique as for the 
actual galaxies (see section [2~Tj) . 

2.1 Stellar mass estimate and volume-limited 
sub-samples 

To calculate the stellar mass for our dataset we use the 
stellar population synthesis results from lBruzual fc Chariot! 
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Figure 1. The distribution of 6dFGS galaxies in log stellar mass 
and redshift. The redshift and stellar mass cuts imposed to create 
the four volume-limited samples (Si to 54 see Table[T} are shown 
by the coloured lines. All galaxies in the upper left quadrant cre- 
ated by the two correspondingly coloured lines are included in the 
volume-limited sub-samples. The plot shows a randomly chosen 
set of 20% of all galaxies. 



(|l995t ) together with a scaled Salpete r initia l mass function 
(IMF) as reported in lBell fc De Jond (|200ll ) 

log 10 (M st ellar/Lj) = -0.57C bj - rp + 0.48 

log 10 (Lj) = (M5 un - Mj)/2.5 (1) 
log 10 (M ste n ar ) = log 10 (A/ stc n ar /ij) +log 10 (Lj), 

with the 2MASS bj-rF colours, Cbj-r F the J-band absolute 
magnitude, Mj an d the J-band a bsolute magnitude of the 
sun, M s / n = 3.70 (|Worthevllr994T ). The biggest uncertainty 
in stellar mass estimates of this type is the choice of the IMF. 
Assuming no trend in IMF with galaxy type, the range of 
IMFs presented in the literature cause uncertainties in the 
absolute normalisation of the stellar M/L-rat io of a factor 
of 2 in the near-infrared (|Bell fc De Jong||200ll ). The 6dFGS 
stellar mass function as well as a comparison of different 
stellar mass estimates is currently in preparation (Jones et 
al. in prep). 

We create four volume-limited sub-samples in redshift 
and stellar mass. This is done by choosing an upper limit in 
redshift (z max ) and then maximising the number of galaxies 
by choosing a lower limit in stellar mass (A^ stc ii ar ), meaning 
that every galaxy above that stellar mass will be detected 
in 6dFGS, if its redshift is below the redshift limit (see Fig- 
ure [TJ. Because of the distribution in bj-rF colour, a clear 
cut in absolute magnitude does not correspond to a clear cut 
in stellar mass and therefore our samples are not perfectly 
volume-limited. We use the absolute magnitude limit which 
corresponds to a chosen redshift limit and derive a stellar 
mass limit, using the bj-v f colour corresponding to the 50% 
height of the bj-r F distribution (bj-r F (50%) = 1.28). We 
create four sub-samples (S1-S4) with upper redshift cuts at 
Zmax = 0.05, 0.06, 0.07 and 0.08 and with the mean log stel- 
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Figure 2. Galaxy density as a function of redshift for the four 
different volume-limited sub-samples S1-S4. All samples follow an 
approximately constant number density indicated by the dashed 
lines and listed in the last column of Table [T] 



lar masses ranging from log 10 (M stc riai-/ft~ 2 Mq) = 10.41 to 
10.73. We also include a low redshift cut-off at z m j n = 0.01. 
Since the stellar mass distributions for these samples over- 
lap, especially for the higher stellar mass samples, the re- 
sults are correlated to some extent (see Figure [I}. All our 
sub-samples are summarised in Table [T] Figure [2] shows the 
galaxy density as a function of redshift for the four sub- 
samples. The roughly constant number density with redshift 
shows that our sub-samples are close to volume-limited. 



From the two-dimensional correlation function, £(r p , 7r) 
we calculate the projected correlation function 



w p (r p ) = 2 / dir£(r p ,ir), 
Jo 



(4) 



where we bin (,(r p ,n) in 30 logarithmic bins from 0.1 to 
IOO/1 -1 Mpc in r p and ty. The upper integration limit in eq.|4] 
was chosen to be 7r max = 50fe _1 Mpc for all sub-samples. We 
tested different values for 7r max , changing it between 20 and 
90/i -1 Mpc without any significant effect to ro or 7 if the er- 
rors on w p (r p ) are adjusted accordingly. The error on scales 
larger than r p = 20 is very large and hence the contribution 
of these scales to the fit is small. 

To derive a covariance matrix for the projected corre- 
lation function we use the method of jack-knife re-sampling 
our galaxy samples. First we divide the dataset into N = 400 
subsets, selected in R.A. and Dec. Each re-sampling step ex- 
cludes one subset before calculating the correlation function. 
The covariance matrix is then given by 



C 1 



(N-l) 
N 



(5) 

where w k (r p ) is the projected correlation function estimate 
at separation rp with the exclusion of subset k and w(r P ) is 
the mean. 

We also note that wide-angle effects can be neglected in 
this an alysis, since we are in tere sted in small-s c ale cl uster- 
ing (see lBeutler et al.l (|201lh and lBeutler et all |2012T ) for a 
detailed investigation of wide-angle effects in 6dFGS). 



3 DATA ANALYSIS 

We measure the separation between galaxies in our sur- 
vey along the line of sight (n) and perpendicular to the 
line of sight (r p ) and count the number of galaxy pairs 
on this two-dimensional grid. We do this for the 6dFGS 
data catalogue, a random catalogue with the same selection 
function, and a combination of data-random pairs. We call 
the pair-separation distributions obtained from this analy- 
sis step DD(r p ,n), RR(r p ,ir) and DR(r p ,n), respectively. 
In the analysis we used 30 random catalogues with the 
same size as the real data catalogue and average DR(r p , it) 
and RR(r p ,n). The random mocks are s a mpled from the 
6dFGS luminosity function (| Jones et alj (|200fJ) and also 
Jones et al, in prep.), and hence they contain the same evo- 
lution of luminosity with redshift that we see in 6dFGS it- 
self. JTJrej^dsliift^spac^^orrelation function is then given by 
the lLandv fc Szalavl (|l993h estimator: 

-1-1. DD< y r v^) ( n r\ 2 I>-R(r p ,7r) ( n r \ 
S[r p ,n) - 1+ RR[rp 7v) 2 RR{rp:7l) {-) - W 

where the ratio n r /rid is given by 

EJV r 

n d E ^ w , ^ 

and the sums go over all random (N r ) and data (Nd) galax- 
ies. Here we employ a completeness weighting, Wi, where 
we weight each galaxy by the inverse sky- and magnitude 
completeness at its area of the sky (Jones et al., in prep.). 



3.1 Fibre proximity limitations 

The design of the 6dF instrumen t does not allow fi bres to 
be placed closer than 5.7 arcmin (| Jones et al.ll2004h . which 
corresponds to a distance of r p « 0.3/i _1 Mpc at redshift 
z — 0.07. This limitation is relaxed in 6dFGS where about 
70% of the survey area has been observed multiple times. 
However, for the remaining 30% we have to expect to miss 
galaxy pairs with a separation smaller than 5.7 arcmin. 

Figure |3] shows the angular correlation function for the 
6dFGS redshift catalogue (w z , blue data points) and the 
target catalogue (w p , solid black line). While the target cat- 
alogue contains 104 785 galaxies, the redshift catalogue of 
galaxies which have a J-band, fej-band and r^-band magni- 
tude contains 76 833. The dashed line indicates the angular 
scale of the fibres. The two angular correlation functions 
agree on scales 6 > 0.1°, with the redshift catalogue falling 
below the target catalogue at lower scales. 

In the lower panel of Figure [3] we show the ratio 
(1 + w p )/(l + w z ). The dashed line shows a spline fit to 
the blue data points, which than can be used to up-weight 
galaxy pairs with small angular separations a nd correct for 
the fibre proximity effect (|Hawkins et al.ll2003f ) . This weight- 
ing is additional to the completeness weighting u>i we intro- 
duced earlier. While the completeness weighting is applied 
to single data- as well as random galaxies, the fibre proxim- 
ity weighting is only applied to data galaxy pairs. Applying 
this weight to the 6dFGS redshift catalogue results in the 
open data points in Figure [3] which now agree very well with 
the target catalogue. 
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Table 1. Summary of the different volume-limited sub-samples, S1-S4, used in this analysis. The effective redshift, z e g and effective 
stellar mass are calculated as the mean of all galaxy pairs contributing to correlation function bins between 0.1 and 40h~ 1 Mpc and the 
error is the standard deviation from the mean. The galaxy density n g is calculated as the number of galaxies N divided by the co-moving 
sample volume, which is calculated from the maximum redshift z max . 

sample w log 10 ( g^ft/A (z) log 10 f ^"""* ) N n 3 [h 3 Mp C - 3 ] 



51 0.05 10.12 0.0369 ± 0.0010 10.410 ± 0.010 24 644 4.536 x 10" 3 

5 2 0.06 10.28 0.0453 ± 0.0011 10.532 ± 0.010 27 999 3.001 x 10" 3 

5 3 0.07 10.42 0.0521 ±0.0012 10.635 ± 0.018 26 584 1.806 X 10" 3 

5 4 0.08 10.55 0.0585 ± 0.0020 10.730 ± 0.013 22 497 1.030 X 10 -3 



r p [h'Mpc] (z eff = 0.07) 
10" 1 1 10 




en 

Figure 3. The angular correlation function of 6dFGS (blue data 
points) as well as the target catalogue (black solid line). The de- 
viation between the two angular correlation functions is caused 
by fibre proximity limitations. The lower panel shows the ratio of 
the two correlation functions, which can be used as a weight to 
correct for the fibre proximity limitations (see section [3.H . The 
open data points show the corrected 6dFGS angular correlation 
function which is in very good agreement with the target cata- 
logue. 



4 POWER LAW FITS 

The projected correlation function can b e related to the 
real-s pace correlation function £(r), using ()Davis fc Peebles! 
Il982ft 

w p (r p )=2 r dyd(r 2 p +y 2 ) 1/2 ] 

roc (6J 

= 2 / rdr£{r)(r 2 - r£) _1/2 . 

J r p 

If the correlation function is assumed to follow a power law, 
£(r) = (r/ro) 7 , with the clustering amplitude ro and the 
power law index 7, this can be written as 

»*> = ^r r (O r (^) /r (i)< m 
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Table 2. Best fitting parameters, ro and 7 for power law fits 
to the projected correlation functions w p (r p ) of the four volume- 
limited 6dFGS sub-samples (S1-S4). The fitting range in all cases 
is 0.1 < r p < 40 h~ 1 Mpc with 24 bins and 2 free parameters. The 
last column shows the reduced \ 2 indicating the goodness of the 
fit. 



sample 


ro [/i _1 Mpc] 


7 


X 2 /cl.o.f. 




Si 


5.14 ± 0.23 


1.849 ±0.025 


23.5/(24 - 2) = 


1.07 


s 2 


5.76 ±0.17 


1.826 ±0.019 


19.5/(24 - 2) = 


0.89 


S3 


5.76 ±0.16 


1.847 ±0.019 


32.6/(24 - 2) = 


1.48 


Si 


6.21 ± 0.17 


1.846 ±0.019 


35.8/(24 - 2) = 


1.63 



with r being the Gamma-function. Using this equation we 
can infer the best-fit power law for £(r) from w p (r p ). 

Table [2] summarises the results of the power law fits to 
the four 6dFGS sub-samples. The best fitting power laws are 
also included in Figure [4] together with the four projected 
correlation functions. We set the fitting range to be 0.1 < 
r p < 40 7i _1 Mpc which includes 24 bins. Although 6dFGS 
has very good statistics at scales smaller than 0.1ft -1 Mpc 
we do not use them for our fits, since at such scales the fibre 
proximity correction becomes more than 30% (see section l3~T1 
and Figure |3J. 

The value of the clustering amplitude, ro, increases from 
5.14 to 6.21/t" 1 Mpc with increasing stellar mass, while there 
doesn't seem to be a clear trend in 7, which varies around 
the value of 7 = 1.84. The reduced x 2 m the last column 
of Table [2] indicates a good fit to the data for the first two 
sub-samples but grows to xVd.o.f. = 1.48 and 1.63 for the 
high stellar mass sub-samples, indicating deviations from a 
power law. The lower panels of Figure [4] show the different 
projected correlation functions divided by the best fitting 
power law (blue data points). Here we can see that the de- 
viations show systematic patterns. Such effects could be re- 
lated to the strong correlations between bins in the correla- 
tion function. Nevertheless, these patterns can be addressed 
with a full HOD analysis, which we will pursue in the next 
section. We will compare the result of our power law fits 
with other studies in a more detailed discussion in section [8] 

We remark at this point, that there is no satisfying the- 
oretical model that predicts a power law behaviour of the 
correlation function, and hence the motivation of such a fit 
can only be empirical. While alternative approaches (e.g. 
the HOD model) rely on assumptions about the clustering 
behaviour of galaxies and dark matter, they are physically 
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Figure 4. The projected correlation function for the four different volume-limited 6dFGS sub-samples (S1-S4). The dashed black lines 
show the best fitting power laws (see Table [2} while the solid black lines show the best fitting HOD models, derived by fitting the data 
between 0.1 ^ r p ^ 40h~ Mpc. The lower panel shows the data and HOD models divided by the best fitting power laws. 



motivated and hence allow us to learn more about galaxy 
clustering than a pure empirical power law fit. 



5 THEORY: HALO OCCUPATION 
DISTRIBUTION 

The Halo Occupation Distribution (HOD) model describes 
the relation between galaxies and mass in terms of the prob- 
ability distribution p(N\M) that a halo of virial mass M 
contains N galaxies of a given type. Knowing how galaxies 



populate dark matter halos, we can use a dark matter corre- 
lation functi on and infer the galax y correlation function. We 
use CAMB (jLewis fc Bridldl2002l ) to derive a model matter 
power spectrum which we turn into a correlation function 
using a Hankel transform 

The underlying cosmological model is fixed to (Qbh 2 , 
n c h 2 , n 3l o-g) = ( 0.0222 7, 0.1116, 0.966, 0.8) as reported 
in iKomatsu et al.l (|201 if ) (we set a$ = 0.9 for one special 



© 0000 RAS, MNRAS 000, 000-000 



6dFGS: Dependence of halo occupation on stellar mass 7 



case). We have to be aware of the fact that the fitted HOD 
parameters depend somewhat on the assumed values of Q m 
and as and hence the absolute values of the HOD param- 
eters could be biased, if the assumed cosmology is wrong. 
However in this study we focus on the relative HOD param- 
eters for different stellar mass selected sub-samples, which 
is fairly robust against such uncertainties. 

Here we employ an analytic HOD meth o dolog y 



2007) 



2010) 



that is similar to that of IZheng, Coil fc Zehavil 
iBlake, Collister fc Lahavl IMPOST ) and IZehavi etal 
We ut ilise analytic appro ximations for the halo mass func 
tion dTinker et al] 120081 s ) . the biased clustering of ha- 



los jTinker et all (20051) and the nonlinear dark matter 
power spectrum (|Smith et al.l 120031 ) . The profile of dark 
mat ter within halos is well descri bed by the NFW pro- 
file (|Navarro. Frenk fc White] ll996T) parameterize d by the 
the concentration-mass relation ( Duffy et al.ll2008l ) 



c(M, 



6.71 



M 



pivot 



1 + *) 



(0) 



with Afpi vot = 2 x 10 12 ft~ 1 M w . We tr i ed rep lacing eq.[5]with 
the form suggested bv lBullock et al.l (1 19991 ). and found that 
the best-fitting HOD parameters changed by much less than 
the statistical errors. 



5.1 HOD framework and formalism 

In the HOD parametrisation it is common to sepa- 
rate the clustering contributions from the most mas- 
sive galaxies, which are assumed to sit in the halo 
centre, from satellite galaxies. This picture of how 



simulations fe.g. iBerlind et al.1 


20031; Simha et al. 


2009) 


and semi-analytic models (e. 


z. IWhite & Frenk 


1991; 


Kauffmann, White & Guiderdoni 


Il993l; Croton et al.l 


200£; 


Bower et al. 20061). The mean central and satellite number 



densit y of galaxies that populate dark matter halos of mass 



densit y ot galaxies that f 
M is jzheng et al.ll2005h 



(Nc(M)) = 
(Ns(M)} = 



if M < M min 

1 if M > M min 

M 
Mi 



(10) 



where M m i n is the minimum dark matter halo mass which 
can host a central galaxy, Mi corresponds to the mass of ha- 
los that contain, on average, one additional satellite galaxy 
({N s (Mi)) — 1) and a sets the rate at which halos accumu- 
late satellites when growing in mass. In very massive halos 
the number of satellites is proportional to halo mass M to 
the power of a. The total HOD number is given by 



(N t (M)) = (iV c (M)> [1 + (N a (M))} 



(11) 



so that a dark matter halo can only host a satellite galaxy if 
it contains already a central galaxy. In our model we assume 
a step like transition from (N C (M)) = to (JV C (M)) = 1. In 
reality this is more like ly to be a gradu al transition with 
a certain width <ri og M dMore et al.l 120091 ') . To account for 
this, other studies (e.g. lZehavi et af]|2010h modify the HOD 
parametrisation of the central halo term to 



(MM)} = - 



1 + erf 



log 10 (M) -log 10 (M g 

0"l og M 



(12) 



which turns into eq.[T0]for the case o"i og M = 0. For our data 
we found that the reduction in x 2 obtained from fits in- 
cluding cicgM as an additional parameter is not big enough 
to justify this parameterisation. For our largest sub-sample 
(S2) we found Ax 2 = -0.78, while a new parameter would 
be justified w hen Ay 2 < — 2 follo wing the Akaike informa- 
tion criterion (jAkaike et al.lll974l ). 

There are also higher parameter models for the satellite 
fraction such as 



(N S (M)) 



M - A/( 
Mi 



(13) 



where Mo is the minimum halo mass at which satellites can 
exist. Again we tested this model and found that the reduc- 
tion in x 2 does not justify this additional parameter. Hence 
we chose the three parameter model (Mi, a, M m i n ) of ca ll (IE 
Within the halo model we can account separately for 
the clustering amplitude of galaxies which sit in the same 
dark matter halo (one halo term) and galaxies which sit in 
different dark matter halos (two halo term). At small scales 
the clustering will be dominated by the one halo term and 
at large scales it will be dominated by the two halo term. 
For the correlation function this can be written as 



f(r) = [l+£ifc(r)]+£ ah (r). 



(14) 



where £ih(r) and ^2h(r) represent the one halo and two halo 
terms respectively. 



5.2 The 1-halo term, $ lh (r) 

We separate the 1-halo term into contributions from central- 
satellite galaxy pairs and satellite-satellite galaxy pairs. The 
central-satellite contribution is given by 



dM ^l Nc (M)N e (M)P^ M) 



g JM vir (r) 



dM 



M 



(15) 

where n g is the galaxy number density and p(r, M) is the 
halo density profile. The lower limit for the integral is the 
virial mass M V ir(r) corresponding to the virial separation 

The satellite-satellite contribution is usually given by a con- 
volution of the halo density profile with the halo mass func- 
tion. Here we calculate this term in k-space since a convo- 
lution then turns into a simple multiplication 

P irW = 4f f~ ' dM^^N c {M)N 2 s {M)\{k,Mf, 

(16) 

where X(k, M) is the normalised Fourier transform of the 
halo density profile p(r,M). The real-space expression of 
eq. (JTeJ) is 



1 



, 2 sin(fcr) 



(17) 



which can then be combined with the central-satellite con- 
tribution to obtain the 1-halo term 



M n 



is fixed by the number density 



(18) 
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5.3 The 2-halo term, &h (r) 

The 2-halo term, ^2h(f), can be calculated from the dark 
matter correlation function since on sufficiently large scales 
the galaxy and matter correlation function are related by a 
constant bias parameter. We calculate the two halo term in 
Fourier space as 



P2h(k,r) =P m (k)x 

•AflimM 



dM \{k,M) 



dM 



n' g (r) 



(19) 



where P m {k) is the non-l inear model power spectrum from 
CAMB including halofit ([Smith et alJl2003l ) and M Um (r) is 
the halo mass limit for which we can find galaxy pairs with 
a separation larger than r v i r . To fix the value of Mij m (r) 
w e use the number de nsity matching model in the appendix 
of iTinker et ail (|2005l ). b h {M,r) is the scale dependent halo 
bias at separation r for which we assume the following 
model (|Tinker et al.ll2005h 



bl(M,r) = bW 1 + 1 - 17Uir)] 



[l + 0.69£m(V 



(20) 



where f m (r) is the non-linear matt er correlation funct ion 
and b(M) being the bias function of ITinker et all (120051 ). 



5.4 Derived quantities 

Since the HOD model is directly based on a description of 
dark matter clustering and its relation to galaxy clustering, 
an HOD model can tell us much more about a galaxy pop- 
ulation, than just the two parameters Mi and a. 

For \(k, M) 1 (which corresponds to large separa- 
tions r in real-space) the two halo term simplifies to 

P2h(k,r)^b 2 cB P m (k), (21) 

where the e ffective bias is the galaxy number weighted halo 
bias factor (|Tinker et alj|2005h 

' dM^^±b h (M)N t (M). (22) 

We can also ask what is the average group dark matter halo 
mass for a specific set of galaxies (often called host-halo 
mass). Such a quantity can be obtained as 



1 

M nff = — 



' dM^MlMN t (M), 



(23) 



which represents a weighted sum over the halo mass function 
with the HOD number as a weight. 

The averaged ratio of satellite galaxies to the total num- 
ber of galaxies is given by 

= J-dM^lN 3 (M) 

13 /" dM^-N c (M) [1 + N s (M)] 

and the central galaxy ratio is given by f c = 1 — / s . 



6 HOD PARAMETER FITS 

To compare the HOD models to the data we use the same 
fitting range as for the power law fits earlier (0.1 < r v < 



40 /i" 1 Mpc). The free parameters of the fit are Mi and a 
as described in the last section and hence we have the same 
number of free parameters as for the power law fits. 

All the fitting results are summarised in Table [3] and 
Figure [4] The reduced \ 2 m the last column of Table [3] in- 
dicates a good fit to the data in all cases and the ratio of 
data to best fitting power law in the lower panels of Fig- 
ure [4] shows that the HOD model reproduces the double 
peak structure present in the data. We also note that the 
reduced \ 2 in case of the HOD fits is uniformly lower than 
for the power law fits indicating a better fit to the data for 
all sub-samples. 

Figure [S] shows the 2D probability distributions in Mi 
and a for the four different volume-limited sub-samples. 
There is a strong trend of increasing Mi with stellar mass 
and a weaker but still significant trend of increasing a with 
increasing stellar mass. This indicates that dark matter 
halos that host a central galaxy with higher stellar mass 
have their first satellite on average at a larger dark mat- 
ter halo mass. However the number of satellites increases 
more steeply with halo mass for samples with higher stellar 
mass. This trend indicates that halos of higher mass have 
greater relative efficiency at producing multiple satellites. 
Si milar trends were fou nd in sub-samples of SDSS galaxies 
by IZehavi et al.l l2005bl . As a test for the sensitivity of the 
fitting results to the upper fitting limit we performed a fit to 
sub-sample 5*2 with the fitting range 0.1 < r v < 20h~ 1 Mpc. 
All parameters agree within la with the fit to the larger fit- 
ting range. We called this fit S'? and included it in Table [3] 

Table [3] also includes derived parameters like the min- 
imum dark matter halo mass, M m i n , the satellite fraction, 
f s , the effective dark matter halo mass, M e ff and the effec- 
tive galaxy bias, b e g. The errors on these parameters are 
calculated as the 68%-confidence level of their ID probabil- 
ity distribution. The effective dark matter halo mass seems 
to be almost constant for all samples while the minimum 
dark matter halo mass M m j n increases with stellar mass. The 
satellite fraction / s decreases with increasing stellar mass, 
indicating that galaxies with high stellar mass have a higher 
probability to be central galaxies. The increasing effective 
galaxy bias indicates that galaxies with higher stellar mass 
are more strongly clustered and hence reside in high density 
regions of the Universe. The increasing effective galaxy bias 
is well described by the following form: 



6(M st ellar) = (1-05 + M^llar/M, ) X (0.8/(7g) 



(25) 



with M* = 1.18 x lO 11 /! 2 M Q . We compare this func- 
tion to the measurements in Figure [6] The absolute stel- 
lar masses are subject to significant uncertainties and dif- 
ferent methods to derive stellar masses can come to very 
different conclusions. However, most stellar mass estimates 
are related by a simple constant offset, and relation 1251 can 
be scaled accordingly. The equation above is only valid for 
the stellar mass range probed in this analysis (M sto iiar = 
2.6 - 5.4 x 10 10 h- 2 M Q ) since the underlying dynamics are 
most likely not captured in eg. 1251 

We will discuss the implications of all these results in 
the next sections. First we will derive 6dFGS mock samples 
using different semi-analytic models. We will then compare 
the predictions from these semi-analytic models to our data 
followed by a comparison to other studies. 



© 0000 RAS, MNRAS 000, 000-000 



6dFGS: Dependence of halo occupation on stellar mass 9 



Table 3. Summary of the best fitting parameters for HOD fits to the four 6dFGS sub-samples (S1-S4) with the fitting range 0.1 < r p < 
40/i — 1 Mpc. The last column shows the reduced x 2 derived from a fit to 24 bins with 2 free parameters and indicates good fits for all 
sub-samples. Error bars on the HOD parameters correspond to lcr, derived from the marginalised distributions. The sample S 2 uses a 
different fitting range of 0.1 < r p < 20fi~ 1 Mpc for the largest of our sub-samples, S2. Furthermore we include a special fit to S2 labeled 
Sj 8 "' 9 , where we change our standard assumption of erg = 0.8 to as = 0.9. The satellite fraction / s . the effective dark matter halo 
mass M e ff and the effective galaxy bias 6 e ff are derived parameters (see section [5.4H . The bias depends on our initial assumption of erg 
and hence this parameter should be treated as b e g X (crg/0.8). 



sample log 10 (■*&—) a log 10 ( , ) fs log 10 ( h -u f r ) &e£E X 2 /d.o.f 



51 13.396 ±0.017 1.214 ±0.031 12.0478 ± 0.0049 0.2106 ± 0.0078 13.501 ±0.015 1.2704 ± 0.0087 13.5/(24 - 2) = 0.61 

5 2 13.568 ±0.013 1.270 ± 0.027 12.2293 ± 0.0043 0.1879 ± 0.0062 13.532 ± 0.012 1.3443 ± 0.0063 16.6/(24 - 2) = 0.75 

53 13.788 ±0.012 1.280 ±0.029 12.4440 ± 0.0039 0.1578 ± 0.0061 13.546 ± 0.012 1.4002 ± 0.0062 21.4/(24 - 2) = 0.97 

5 4 14.022 ±0.011 1.396 ±0.033 12.6753 ± 0.0032 0.1243 ± 0.0058 13.617 ± 0.012 1.5170 ± 0.0077 24.4/(24 - 2) = 1.11 



££8-0.9 i3. 6 05 ±0.013 


1.193 ± 0.023 


12.2382 ±0.0048 


0.2034 ±0.0077 


13.649 ±0.014 


1.2151 ±0.0065 


19.6/(24 


-2) 


= 0.89 


S' 2 13.568 ±0.013 


1.282 ± 0.027 


12.2291 ± 0.0042 


0.1890 ±0.0074 


13.542 ±0.010 


1.3492 ±0.0065 


19.1/(21 


-2) 


= 1.01 



7 SEMI-ANALYTIC MOCK CATALOGUES 

To compare our result with theory, we derive 6dFGS 
moc k catalogues from two different semi -analytical mod- 
els i|Croton et al.1 12006J ; iBower et all 120061). both based on 
the Millennium Simulation ( Springel et all 120051 ) publicly 
available through the Millennium Simulation databasq 2 ]. 
Semi-analytic models are based on an underlying N-body 
simulation together with theoretically and observationally 
motivated descriptions of gas cooling, star formation and 
feedback processes. 

The Millennium Simulation is a dark matter only N- 
body simulation which traces the hierarchical evolution of 
2160 3 particles in a periodic box of 500 3 /i~ 3 Mpc 3 from 
redshift z — 127 to z = 0. The un derlying cosmologi - 
cal model follows WMAP1 cosmology (|SpergeI et all 120031 ) 
given by a matter density of f2 m = fldm + = 0.25, 
a cosmological constant of Qa = 0.75, a Hubble con- 
stant of Ho = 75km/s/Mpc, a spectral index of n s = 1 
and a r.m.s. of matter fluctuations in 8/i _1 Mpc spheres of 
<rg = 0.9. The individual particle mass of the simulation 
is 8.6 x 10 8 h' 1 Mq and halos and sub-halos are identified 
from the spatial distribution of dark matter particles using a 
standard friends-of-friends al gorithm and the SUBFIND al- 
gorithm (|Springel et alll200ll ). All sub-halos are then linked 
together to construct the halo merger trees which represent 
the basic input of the semi-analytic models. 

Here we are using the 2 = output of the Millennium 
Simulation. To ensure that the stellar masses are calculated 
in a consistent and comparable way, we implement the fol- 
lowing procedure to derive the 6dFGS mock catalogues: 

(i) We apply the 6dFGS Jf-band apparent magnitude lim- 
its of 8.85 < K < 12.75 to the full 500/i" 3 Mpc 3 simulation 
box. We have to use the Tv-band instead of the J-band, 
which is actually used in this analysis, because none of the 
semi-analytic models provide J-band magnitudes. However, 
the 6dFGS J-band and 7v"-band samples have significant 
overlap, meaning that almost all galaxies which have a K- 



2 http:/ /gavo. mpa-garching.mpg.de/MyMillennium/ 
http: / /galaxy-catalogue. dur.ac.uk:8080/Millcnnium/ 



band magnitude also have a J-band magnitude. We also 
account for sky- and magnitude incompleteness. 

(ii) We re-calculate stellar masses for each galaxy using 
the technique described in section | 2. II but with t he corre- 
sponding K-h&nA relations (|Bell fc De Jong||200ll ). instead 
of J-band. The re-calculation of the stellar masses ensures 
that potential disagreement with our measurement is not 
caused by a different technique of deriving stellar masses or 
a different assumptions about the IMF. 

(iii) We apply the same redshift and stellar mass limits 
to the semi-analytic catalogues, which we used to produce 
the four volume-limited samples in the 6dFGS dataset (see 
Table EE). 

We found that all mock 6dFGS catalogues derived from 
these semi-analytic models contain fewer galaxies than the 
data sample (by about 40%). This could be related to the 
slightly different cosmology used in these simulations, which 
should have its largest impact on large clusters, which are 
sampled in 6dFGS. We are not attempting to correct for 
such differences in the cosmological model. The aim of this 
part of our analysis is to test the current predictive power 
of semi-analytic models. 

Semi-analytic models are often gr ouped into "Durham 
models" and "Munich models". The IBower et all (|2006l ) 
mo del belongs in the g roup of "Durham models" while 
the ICroton et al. I (|2006l ) model belongs in the group of 
"Munich models". 



7.1 Durham models 

In the Durham mode ls, merger trees are produced follow- 
ing jHeUxie^L|_j200^)which_are independent of those gen- 
erated by Springel et all (|2005l ). When the satellite galaxy 
falls below a certain distance to the central galaxy given by 



r c + r s 



where r c and r s are the half mass radii 



of the central and satellite galaxy, respectively, the satel- 
lite and central galaxy are treated as one. The largest of 
the galaxies contained within this new combi ned dark mat- 
ter h alo is assumed to be the central galaxy (|Benson et al.l 
120021) . whilst all other galaxies within the halo are satellites. 
The dynamical friction and tidal stripping which are present 
in such a system are modelled analytically. These models are 
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Figure 5. Two dimensional probability distribution for 
log 10 (Mi) and a for the HOD fits to the four 6dFGS volume- 
limited sub-sample (S1-S4). We include two fits for the largest 
sample 52, where for S^ 8 we change our standard assump- 
tion of ag = 0.8 to erg = 0.9. The parameters derived from the 
fits are summarised in Table [3] The best fitting values arc marker 
with black crosses. 



• 6dFGS (this analysis) 




I 1.05 + M m||lt /1.18x10"hX 









10.3 10.35 10.4 10.45 10.5 10.55 10.6 10.65 10.7 10.75 10.8 
log (M /h 2 WU 

M 10 v stellar 



Figure 6. The effective galaxy bias (see eq. I22H as a function of 
log stellar mass for 6dFGS (blue data points). The increase in 
galaxy bias with stellar mass can be described by eq. I25I which 
corresponds to the black line. 



based on NFW density profiles for the central halo as well 
as the satellite halo, while galaxies are modelled as a disc 
plus spheroid. 

7.2 Munich models 

T he Munich models ar e based on the original merger trees 
by ISpringel et all (|2005t ) . One of the key differences between 
these merger trees and the ones used in the Durham models 
is that the Munich models explicitly follow dark matter halos 
even after they are accreted onto larger systems, allowing 
the dynamics of satellite galaxies residing in the in-falling 
halos to be followed until the dark matter substructure is 
destroyed. The galaxy is than assigned to the most bound 
particle of the sub-halo at the last time the sub-halo could 
be identified. 

In the Munich models, a two-mode formalism is adopted 
for active galactic nuclei (AGNs), wherein a high-energy, or 
"quasar" mode occurs subsequent to mergers, and a con- 
stant low-energy "radio" mode suppresses cooling flows due 
to th e interaction betw een the gas and the central black 
hole (|Croton et al.ll200fj) . In the quasar model, accretion of 
gas onto the black hole peaks at z ~ 3, while the radio mode 
reaches a plateau at z ~ 2. AGN feedback is assumed to be 
efficient only in massive halos, with supernova feedback be- 
ing more dominant in l ower-mass halos. 

The galaxies in the lCroton et al] (|2006l ) model can be of 
three different types: central galaxies (type 0), satellites of 
type 1 and satellites of type 2. Satellites of type 1 are associ- 
ated with dark matter substructures, which usually refers to 
recently merged halos. Satellites of type 2 are instead galax- 
ies whose dark matter halo has completely merged with a 
bigger halo and are not associated with a substructure. We 
treat both type 2 and type 1 as satellite galaxies. 

7.3 Testing semi-analytic models 

In their original paper iBower et al.l (|2006l ) compare the 
K-band luminosity function, galaxy stellar mass function, 
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and cosmic star formation rate with high-redshift obser- 
vations. They find that their model matches the observed 
mass and luminosity functi ons reasonably well up to z « f . 
iKitzbichler fc White! i|2007l ) compare the magnitude counts 
in the bj, r\p, / and K -bands, redshift distributions for K- 
band selected samples, bj- and K-band lu minosity func- 
tions, and galaxy stellar mass functio n from Croton et al 



2006 ) and the very similar model by IDe Lucia fc Blaizot 



20071 ) with high redshift measurements. They find that the 
agreement of these models with high-redshift observations 
is slightly worse than that found for the Durham models. 
In particular, they find that the Munich models tend to sys- 
tematically overestimate the abundance of relatively massive 
gala xies at high r e dshift . 

ISnaith et al.l (l201lf) compared f o ur d if ferent semi- 
analytic models (IDe Lucia fc Blaizo ij 120071 ; | Bower et al.l 
12009 ; IBertone et alj|2007l ; IFont et all 120081 '). with observa- 
tions and found that all models show a shallower, wider 
magnitude gap, between the brightest group galaxy and the 
seco nd brightest, com p ared t o observations. 

Ide la Torre et all J2010TI compared mea s ureme nts of 
VVDS with the model bv lDe Lucia fc BlaizotI (|2007f ). They 
found that the model reproduced the galaxy clustering at 
2 > 0.8 as well as the magnitude counts in most bands. How- 
ever the model failed in reproducing the clustering strength 
of red galaxies and the bj — I colour distribution. The 
model tends to produce too many relatively bright red satel- 
lites galaxies, a fact that has been reported in other stud- 
ies as well (over-quenching proble m: IWeinmann et al1l2006l : 
iKimm et alfaod ; iLiu et al)|201Ch . 



8 DISCUSSION 

8.1 Effective halo mass and satellite fraction 

As is evident from Table [3] the 6dFGS galaxies sit in mas- 
sive central dark matter halos and most of our galaxies are 
central galaxies in these halos with only a small fraction 
being satellite galaxies. Because of the way 6dFGS galax- 
ies are selected, the large majority are red elliptical galaxies 
and hence our findings agree very well with previous studies, 
which also found that such galaxies are strongly clustered 
and therefore must reside in high density regions. 

First we will discuss the effective dark matter halo mass, 
which represents the effective group or cluster mass for the 
distribution of galaxies (not the mass of the individual ha- 
los which host the galaxies). This quantity appears almost 
constant for our four different sub-samples. The increasing 
central galaxy mass (M m i n ) with stellar mass is offset by 
the decreasing satellite fraction, resulting in a fairly con- 
stant M c ff. While the HOD model allows us to derive only 
averaged parameters for each sample, semi-analytic models 
directly connect dark matter halo masses with single galax- 
ies. Figure [7] shows the distribution of central and satellite 
galaxies as a fu nction of friends-of- f riends (fof ) halo mass de- 
rived from the iBower et all (|2006T ) semi-analytic catalogue 
together with the 6dFGS selection criteria. The catalogue 
contains fewer satellites than central, but the satellites sit 
in very massive dark matter halos and hence have a sig- 
nificant impact on M e ff. This plot shows that most 6dFGS 
galaxies sit in fO 11 — f0 12 ft _1 A/0 halos, while the satellites 
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Figure 7. This plot 
mass login(Mf n f ) of the 



shows the friends-of-friends halo 
satellite and central galaxies from 



the IBower et al. [ 120061) semi-analytic catalogue with the 6dFGS 
selection. The bimodal distribution is caused by a dominant frac- 
tion of halos, in which the 6dFGS galaxies are the central galaxy 
and a smaller fraction with larger halo mass, in which the 6dFGS 
galaxies are satellite galaxies. 



Table 4. Summary of parameters derived from different semi- 
analytic models. Wc impose the 6dFGS _R"-band apparent mag- 
nitude selection (J-band is not available for these semi-analytic 
models) as well as correct for incompleteness. The stellar masses 
are re-calculated using the technique described in section \2. II To 
calculate the effective dark matter halo mass we used the friends- 
of-friends halo mass. The stellar mass and the effective halo mass 
arc calculated as the mean of the distribution. 



sample 



fs log! 



' (M a , 



log! 



f Croton et al. 20061 



SI 


0.138 


10.38 


13.400 


S2 


0.099 


10.48 


13.389 


S3 


0.085 
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13.507 


SI 


0.090 


10.67 


14.013 


(Bower et al. 2006) 


SI 


0.28 


10.36 


13.695 


S2 


0.25 


10.47 


13.661 


S3 


0.22 


10.58 


13.674 


S4 


0.18 


10.68 


13.777 



sit in very massive groups and clusters of up to 10 h~ Mq. 
While the median would be around 10 12 /i -1 Mq, the effec- 
tive mass M e g from the HOD model is the mean of this 
distribution, which is pushed to very large values by the 
satellite fraction. 

We calculated the effective halo mass M e e — {Mf f) 
for the four volume-limited sub-samples derived from the 
two semi-analytic models and summarise these results in Ta- 
ble [4] While the effective h alo mass ap pears constant i n case 
of the IBower et ail l|200rj ) model, the ICroton et all (|2006t ) 
model shows an increase with increasing stellar mass. These 
different behaviours for the two different semi-analytic mod- 
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Figure 8. Satellite fraction of the 6dFGS stellar mass 
sub-samples (blue) and the corresponding semi-analytic sub- 
samples based on the Millennium simulation llCroton et al.ll2006t 
iBower et alj|2006h . The red data point represents a fit to sample 
S2 which assumes ag = 0.9, while the blue data points assume 
as = 0.8. 



els are most likely connected to the different trend in the 
satelli te fraction shown in Figure [8] While the lBower et al.l 
(|2006l ) model shows a constant decrease in the satellite frac- 
tion, the ICroton et al.l ([2006) model reaches a constant at 
large stellar mass which causes the effective halo mass to 
rise. A lower satellite fraction at a fixed stellar mass means 
that at fixed halo mass the s atellite galax i es ar e less mas- 
sive |Weinmann et alj 120061 ; iKimm et all 120091 : iLiu et all 
l2010h . 

F rom Figure [5]we can se e that both the ICroton et all 
(|200rJ ) and the IBower et all (|200rJ ) model give slightly 
smaller stellar masses than measured in 6dFGS. This in- 
dicates that the semi-analytic catalogues contain too few 
bright galaxies at low redshift. 

Figure [8] also includes the fit to sample S2 which as- 
sumes as — 0.9 (red data point), which agrees with the as- 
sumptions of the Millennium Simulation, while the blue data 
points assume as = 0.8. A larger erg causes a larger satellite 
fraction which than causes a larger effective hal o mass. This 
effect seem to bring our results closer to the IBower et al.l 
|2006r i model. 



8.2 The Mi - M min scaling relation 

The HO D analysis of early data from SDSS bv lZehavi et al.l 
(2005b) showed the relation, Mi « 23M m m, between the 
mass of halos that on average host one additional satellite 
galaxy, Mi and the minimum dark matter halo mass to host 
a central galaxy, M m i n . This h as been confirmed in subse- 
quent studies (Mi sa 18 M min inlZheng, Coil fc Zehavill2007l 
and Mi « 17M„,, JZehavi et alj|20lj ). This relation implies 



that on average a halo hosting two galaxies of the type stud- 
ied in their analysis has a mass ~ 23 times the mass of a 
hal o hosting only one g alaxy of the same type. 

IZehavi et al.l (|2010T l also found that this scaling factor is 
somewhat smaller at the high luminosity end, corresponding 
to massive halos that host rich groups or clusters. This latter 



trend likely reflects the relatively late formation of these 
massive halos, which leaves less time for satellites to merge 
onto central galaxies and thus lowers the satellite threshold 
Mi. 

Theo retical studies of HO Ds in dark matter sim- 
ulations (|Kravtsov et al.l 12004 ) and those predicted by 
Smoothed-particle hydro dynamics (SPH) a nd semi-analytic 
galaxy formation models (|Zheng et al.ll200"5l ) reveal a similar 
relation with a scaling factor of ss 20. The large gap between 
Mi and M m i n arises because in the low occupation regime, 
a more massive halo tends to host a more massive central 
galax y, rather than multiple smaller galaxies |Berlind et al.l 

l2003h. 

" TAbbas et~all (|2010h did a similar study using data from 
the /-band selected VIMOS-VLT Deep Survey (VVDS) and 
found the ratio Mi/M m i n to be ~ 40 — 50. The galaxies in 
this study are at much higher redshift (z « 0.83) compared 
to the SDSS galaxies. This result means that in order to 
begin hosting satellite galaxies, halos sampled by the VVDS 
survey need to accrete a larger amount of mass compared to 
SDSS halos. 

Wake et al.l (|201ll ) also studied high redshift galaxies 
(1.1 < z < 1.9) in the NEWFIRM Medium Band Survey 
(NMBS) and found the rati o Mi/M min to be « 4-10. This, 
together with the result by I Abbas et~afl i]201Ch shows, that 
this relation is not as fundamental as originally thought, 
but strongly dependent on the type of halos probed in each 
ana lysis. 

iLeauthaud et all (|2012j ) studied galaxy clustering in the 
COSMOS survey using threshol d stellar mass samples. Us- 
ing the equations discussed in |Leauthaud et af]l201ll ) we 
can derive Mi and M m i n from their HOD parametrisation, 
whi ch is included in F igure [9] 

iMatsuoka et all (|201lh analysed ~ 60 000 massive 
(log lo (M st eiiar//i" 2 M ) > 10.7) galaxies from the UKIRT 
Infrared Deep Sky Survey (UKIDSS) and the SDSS II 
Supernova Survey. This analysis shows a very different 
clustering amplitude depending on whether the observed 
and theoretical number densities are matched up or not. 
This dependency makes a comparison of our results with 
their derived parameters very difficult. Nevertheless, the 
ratio Mi /M m i n in their analysis does not depend signifi- 
cantly on their initial assumptions and hence we included 
their results in Figure |^|. Their results indi cate a lower 
Mi /M m j n ratio at large Mi consistent w ith IZehavi et al.l 
J2010l ) and lBlake. Collister fc Lahavl (|2008l ). 

In our analysis we found a scaling relation of Mi w 
22 M m in (the exact values are 22.29 ± 0.39,21.81 ± 
0.30, 22.08 ± 0.28 and 22.22 ± 0.25 for S1-84, respectively). 
In Figure [9] we compare our results with other stellar mass 
selected samples (coloured data points) and with luminos- 
ity threshold selected samples (black data poin ts) . The ratio 
Mi /Mm in found in 6dFGS is in agree ment with | Zehavi et al.l 
(|2005tJ) and sligh t ly larg er t han IZehavi et al.l (|2010T ). 
While IZehavi et all (l2005bh and IZehavi et al l (|2010f l study 
the same type of galaxies. IZehavi et al. ( 2005bl ) uses an HOD 
param eterisation very similar to ours, while IZehavi et al.l 
(|2010l) uses a parameterisation based on 5 free parameters, 



3 We use their results, in which they match the observed and 
theoretical number densities, since this agrees with our method. 
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Fi gure 9. The relatio n between M± and M m j n for the four different volum e-limited 6dFGS sampl es (blue data point s), compared 
to IZehavi et al.| l|2005bl) (blac k solid squares) , iBlake. C ollister &: Lah avl j2008T l (black open triangle s). IZeha vi et al.| | 2010l ) (black open 
squares). lAbbas et alj d2010l) (black solid triangles). iMatsuoka^t al.l ll20 1ll) ( black open points), IWake et"lrtT "i l201ll) (magenta open 
triangles) and iLeeiuthaud et al.l l|2012l ) (magenta solid triangles). IZehavi et alj l|2005bl ) don't report errors on their parameters and the 
6dFGS errors are smaller than the data points. All studies which investigate HOD as a function of stellar mass are coloured (including 
our analysis), while luminosity defined studies are in black. 



suggesting that the differences might be related to the pa- 
rameterisation. 



8.3 Comparison to iMandelbaum et all (|2006h 

An alternative method for probing the connection between 
stellar mass and halo mass is galaxy-g alaxy weak lensing. 
Here we will compare our findings to IMandelbaum et al.l 
(|2006l ) , who used weak lensing of ~ 350 000 galaxies from 
SDSS and looked at the dependence of the amplitude of 
the lensing signature as a function of galaxy type and stel- 
lar mass. They used stellar masses der ived from the z- 
band magnitude and the ratio M/L z fromlKauffmann et al.l 
(|2003r ) with the assumption of a lKroupal (|200ll ) IMF. 

The HOD model employed in there analysis is given by 



N t (M) = N C (M)+N S (M) 



(26) 



with the central galaxy number given by eq. [TO] while the 
satellite galaxies number is modelled by a step-like function 



N S (M) 




kM 

kM 2 



if M > 3M min 

if M min < M < 3M min (27) 
if M < M min . 



The normalisation constant k can be determined by match- 
ing the measured satellite fractions /™ a and is given 



in appendix [X] Here we have adjusted the nomencla- 
ture to the one used in o ur analysis. We also note 
that IMandelbaum et all (|2006l ) assumed as = 0.9, while we 
assumed as — 0.8 for most of our fits. 

From the equation above we can derive the effective 
halo mass using eq. [231 together with the dark matter halo 
mass function. 

We compare the repor ted satellite fraction, /" 



of IMandelbaum et al.1 (|2006i ) and the derived effective 
halo mass, M c g, with our results in Figure [TO] and 
Fi gure 1 1 1 1 Figure 1101 show s that the galaxies sampled 
by Mandelbaum et al. (|2006T ) have a very similar effective 
halo mass compared to 6dFGS galaxies. The 6dFGS mea- 
surements follow the trend of early-type galaxies in these 
measurements, while the late types have a slightly smaller 
eff ective halo mass. However, the errors in case of the results 
by IMandelbaum et al.l (|2006l ) don't allow one to distinguish 
between the early and late type galaxies. Since effective halo 
mass and satellite fraction are strongly linked, we also com- 
pare the satellite fraction of lMandelbaum et al.l |2006i ) with 
our results in Figure 1111 Both figures also include the fit 
to sample S2 where we assumed as = 0.9, instead of the 
st andard as = 0.8, since th is agrees with the assumption 
in IMandelbaum et al.l (|2006l ) . A larger as increases the ef- 
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Figure 10. The effective dark matter halo mass of the 
four 6dFGS sub-samples (blue data points) compared to early 
type (black triangles) and lat e type (black squares) galaxies 
from IMandclbau m et al In red wc also include the fit 

to sample S2 which assumes ag = 0.9 i n agreement with the 
assumptions of IMandclba um et al.l |200li I, while the blue data 
points assume erg = 0.8. 
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Figure 11. The satellite fraction of the four 6dFGS sub-samples 
(blue data points) compared to early type (black triangles) an d 
late type (black squares) galaxies from Mandelb aum et all (2006). 
In red we also include the fit to sample 52 which assumes erg = 0. 9 
in agreement with the assumptions of iMandclb aum et al .1 (2006), 
while the blue data points assume <rg = 0.8. 



fective halo mass by a small amount and hence does not 
impact this comparison significantly. 

Overall we see ve ry good agreement betwee n our re- 
sults and the results of iMandelbaum et al.l (|200fj| ). which is 
reassuring since the two techniques are subject to different 
systematic uncertainties. We can also emphasise that the 
6dFGS results are extremely precise compared to the lens- 
ing results. 



8.4 Comparison to other studies 

iMeneux et al.l (|2007l ) study galaxy cluster mass dependence 



on stellar mass in the VIMOS-VLT Deep Survey (VVDS) 
at redshift 0.5 < z < 1.2. The stellar masses in this sam- 
ple cover the range log lo (M stellar //i- 2 M ) = 8.7 - 10.7. To 
quantify the clustering, they use power laws fits to the pro- 
jected correlation function w p (r p ) and found an evolution in 
r from 2.76 ± 0.17/T 1 Mpc at log 10 (M st eiiar//i~ 2 M©) > 8.7 
to ro = 4.28 ± 0.45ft" 1 Mpc at log 10 (M ste ii a r//i" 2 M Q ) > 
10.2. The slope changes over the same range from 7 = 
1.6 7 ± 0.08 to 7 = 2 .28 ± 0.28. Compari ng the results 
by IMeneux et al] \200H ) with lLi et ail (|2006l ) who used the 
clustering of 200 000 SDSS galaxies at z = 0.15, showed 
that the evolution of the amplitude and shape of the cor- 
relation function w p (r p ) with redshift is faster for low stel- 
lar mass objects than for high stellar mass objects. At low 
stellar mass, the amplitude of ui p (r p ) increases by a factor 
~ 2 — 3 from high to low redshift, while at the high stellar 
mass range log 10 (M s t e iiar /h~ 2 Mq) = 10.2 — 10.7, the ampli- 
tude at 2 ~ 0.85 and z ~ 0.15 are very similar, within the 
error bars. In 6dFGS we found a significantly larger clus- 
tering amplitude, ranging from r = 5.14 ± 0.23/i _1 Mpc at 
log lo (M stollar //r 2 M ) = 10.41 to r = 6.21 ± O.m^Mpc 
at log 10 (M s tcii ar /' ! i _2 A/Q) = 10.73, while our power law in- 
dex 7 stays constant at 7 ~ 1.84. This shows that the trend 
observed in IMeneux etail (|2007l ) and iLi et al.l |2006j) con- 
tinues to the 6dFGS redshifts, even for high stellar mass 
galaxies. 



9 CONCLUSIONS 

We present in this paper an analysis of the clustering proper- 
ties of four stellar mass selected volume-limited sub-samples 
of galaxies from the 6dF Galaxy Survey. The stellar masses 
are calculated using the J-band ma gnitude and the b j — t f 
colour following the technique by iBell fc De Jond (|200lf ). 
The average log-stellar mass for the different sub-samples 
ranges from log lo (M stellar /ft- 2 M ) = 10.41 - 10.73. Our 
analysis has the following main results: 

• The projected correlation function, w p (r p ), for the two 
low stellar mass sub-samples (Si and 5*2) can be described 
by a power law with an acceptable x 2 , while the two high 
stellar mass sub-samples (S3 and S4) have a reduced x 2 of 
1.48 and 1.63, respectively. We also find patterns in the de- 
viations between the best power law fit and the data, which 
can naturally be explained within the halo model. This is 
in agreement with theoretical studies (| Watson et al.l[201ll ). 
predicting that the disagreement of w p (r p ) with a power law 
fit should grow with clustering amplitude. 

• We used an HOD parameterisation with three free pa- 
rameters (Mi, a, M m i n ), representing the typical halo mass, 
satellite power law index and minimum host halo mass, re- 
spectively. The minimum halo mass is fixed by the galaxy 
number density for all our parameter fits. We performed 
fits to the projected correlation functions of the four differ- 
ent volume-limited sub-samples. We tested alternative HOD 
parameterisations, but found that our data does not jus- 
tify more free parameters. Mi and a show increasing trends 
for increasing stellar mass, with log 10 (Mi//i~ 1 M@) rang- 
ing from 13.4 — 14 and a from 1.21 — 1.4. This means that 
galaxies with larger stellar mass populate larger dark matter 
halos, which accrete satellites faster with increasing mass, 
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compared to dark matter halos populated by galaxies with 
lower stellar mass. 

• From the halo model, we can derive averaged param- 
eters for the four samples such as the satellite fraction, ef- 
fective dark matter halo mass and the effective galaxy bias. 
We found that the satellite fraction decreases with stellar 
mass from 21% at log 10 (M ate n al .//i- 2 M ) = 10.41 to 12% 
at log 10 (M sto n ar //i -2 M Q ) = 10.73. The effective dark mat- 
ter halo mass stays constant at log 10 (M e ff//i -1 Mo) « 13.55 
for all four sub-samples. The effective galaxy bias increases 
with stellar mass indicating that galaxies with higher stellar 
mass reside in denser environments. The increase in the ef- 
fective galaxy bias can be described by (1.05+M stG ii al ./M*) x 
(ag/0.8) with M* = 1.18 x lO n /i" 2 M . 

• We use two sem i-analytic models, b ased on the Mil- 
lenni um Simulation (jCroton et al.l 120061 and iBower et alj 
120061 ) to derive 6dFGS mock surveys. We compare the re- 
sults oFthsse_s^mi-analytic models with our measurements. 
The lCroton et all d2006h model u nder- predicts the satellite 
fraction, while the Bower et al.l (|200rj) model is in better 
agreement with our observations, although it slightly over- 
predicts the satellite fraction. Since the effective dark mat- 
ter halo mass is stron gly linked to the satellite fraction, 
the IBower et al.l (|2006l ) model prediction of M e g is again 
in better agreement with our observations. 6dFGS allows a 
powerful test of semi-analytic models, because of (1) the 
robust stellar mass estimates and (2) the focus on "red 
satellites", which semi-analytic models struggled to repro- 
duce in the past. Our results can be used as a new constraint 
on semi-analytic models in the future. 

• For the four volume-limited samples we find a constant 
scaling relation between Mi and M m i n of M i 22M m j n , 
in agreemen t with studies of SDSS galaxies |Zehavi et al.l 
I2005bl , l20ioh . This indicates that 6dFGS galaxies populate 
dark matter halos in a similar way to SDSS galaxies. How- 
ev er, we see a w i de spe ctrum of this r atio ranging from > 40 
in lAbbas et all (|2010h to < 10 in i|Matsuoka et alj|201ll ; 
IWake et al. 20111 ) indicating that this ratio is not univer- 
sal for all types of galaxies. 

• We compare o ur r esults with the results 
of iMandelbaum et al.l (|2006h from galaxy-galaxy weak 
lensing. We find overall good agreement which represents 
a valuable crosscheck for these two different clustering 
measurements. Although our analysis depends on slightly 
more assumptions and covers a smaller range in stellar 
mass, the 6dFGS results are extremely precise compared to 
the lensing results. 
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APPENDIX A: NORMALISATION CONSTANT 
IN Mandelbaum et al. 2006 

We can determine the normalisation constant k of eg. 1271 by 
matching to the satellite fraction /™ a , which in terms of the 
mass function dn(M) / dM is given by 

r a = S °°dM^lN B (M) 

Jo°° dM± j£r [Nc{M) + N s (AO] 



kH 



with 



SudM^ + kH 



1 f min ,* r dn(M) ,, a 
H =—— — / dM — tt-t-M 2 + 

3Af min J M . dM 



which leads to 



,»,dn(M) „ 



rma r°o j . r dn(M) 
Js JAf„i„ UU1 dM 



(A2) 



(A3) 



1 — T^-TT-- M 

where k has units of 1/Ii~ 1 Mq and is constant for a given 
set of / s ma and M min . 
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